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We show that the Quantum Spin Hall Effect, a state of matter with topological 
properties distinct from conventional insulators, can be realized in HgTe/CdTe 
semiconductor quantum wells. By varying the thickness of the quantum well, 
the electronic state changes from a normal to an "inverted" type at a critical 
thickness d c . We show that this transition is a topological quantum phase tran- 
sition between a conventional insulating phase and a phase exhibiting the QSH 
effect with a single pair of helical edge states. We also discuss the methods for 
experimental detection of the QSH effect. 
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The spin Hall effect (|7]|21I21E1I21) has attracted great attention recently in condensed matter 
physics both for its fundamental scientific importance, and its potentially practical application 
in semiconductor spintronics. In particular, the intrinsic spin Hall effect promises the possibility 
of designing the intrinsic electronic properties of materials so that the effect can be maximized. 
Based on this line of reasoning, it was shown (6) that the intrinsic spin Hall effect can in prin- 
ciple exist in band insulators, where the spin current can flow without dissipation. Motivated 
by this suggestion, the quantum spin Hall (QSH) effect has been proposed independently both 
for graphene ( 7) and semiconductors (0|9]) , where the spin current is carried entirely by the 
helical edge states in two dimensional samples. Time reversal symmetry plays an important 
role in the dynamics of the helical edge states ((7011771172]) . When there are an even number of 
pairs of helical states at each edge, impurity scattering or many-body interactions can open a 
gap at the edge and render the system topologically trivial. However, when there are an odd 
number of pairs of helical states at each edge these effects cannot open a gap unless time rever- 
sal symmetry is spontaneously broken at the edge. The stability of the helical edge states has 
been confirmed in extensive numerical calculations (173) IT?)) . The time reversal property leads 
to the Z 2 classification (ITOl) of the QSH state. States of matter can be classified according to 
their topological properties. For example, the integer quantum Hall effect is characterized by a 
topological integer n (|75l) . which determines the quantized value of the Hall conductance and 
the number of chiral edge states. It is invariant under smooth distortions of the Hamiltonian, as 
long as the energy gap does not collapse. Similarly, the number of helical edge states, defined 
modulo two, of the QSH state is also invariant under topologically smooth distortions of the 
Hamiltonian. Therefore, the QSH state is a topologically distinct new state of matter, in the 
same sense as the charge quantum Hall effect. 

Unfortunately, the initial proposal of the QSH in graphene ( 7) was later shown to be unreal- 
istic di^O, as the gap opened by the spin-orbit interaction turns out to be extremely small, of 
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the order of 1CT 3 meV. There are also no immediate experimental systems available for the pro- 
posals in Ref. ([#][7#]). Here, we present theoretical investigations of the type-Ill semiconductor 
quantum wells, and show that the QSH state should be realized in the "inverted" regime where 
the well thickness d is greater than a certain critical thickness d c . Based on general symmetry 
considerations and standard k ■ p perturbation theory for semiconductors (1791) . we show that 
the electronic states near the T point are described by the relativistic Dirac equation in 2 + 1 
dimensions. At the quantum phase transition at d = d c , the mass term in the Dirac equation 
changes sign, leading to two distinct {7(1) -spin and Z 2 topological numbers on either side of 
the transition. Generally, knowledge of electronic states near one point of the Brillouin Zone 
is insufficient to determine the topology of the entire system, however, it does give robust and 
reliable predictions on the change of topological quantum numbers. The fortunate presence of 
a gap closing transition at the T point in the HgTe/CdTe quantum wells therefore makes our 
theoretical prediction of the QSH state conclusive. 

The potential importance of inverted band gap semiconductors like HgTe for the spin Hall 
effect was pointed out in (019]). The central feature of the type-Ill quantum wells is band 
inversion: the barrier material such as CdTe has a normal band progression, with the T 6 s- 
type band lying above the r§ p-type band, and the well material HgTe having an inverted band 
progression whereby the s-type T 6 band lies below the p-type r 8 band. In both of these materials 
the gap is the smallest near the T point in the Brillouin zone (Fig. 1). In our discussion we 
neglect the bulk split-off T 7 band, as it has negligible effects on the band structure (12011271) . 
Therefore, we shall restrict ourselves to a six band model, and start with the following six basic 
atomic states per unit cell combined into a six component spinor: 

vf = (|r 6 , |r 6 , -\), |r 8 , |>, |r 8 , |r 8 , -|>|r 8 , -|>) . (i) 

In quantum wells grown in the [001] direction the cubic, or spherical symmetry, is broken 
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down to the axial rotation symmetry in the plane. These six bands combine to form the spin 
up and down (±) states of three quantum well subbands: El, HI, LI (|2iT) . The LI subband 
is separated from the other two ([271) . and we neglect it, leaving an effective four band model. 
At the T point with in-plane momentum ki\ = 0, mj is still a good quantum number. At this 
point the \E1, mj) quantum well subband state is formed from the linear combination of the 
|r 6 , mj = ±|) and the |r 8 , mj = ±|) states, while the \H1, mj) quantum well subband state 
is formed from the |r 8 , mj = ±|) states. Away from the T point, the El and the HI states can 
mix. As the |r 6 , mj — ±|) state has even parity, while the |r 8 , mj — ±|) state has odd parity 
under two dimensional spatial reflection, the coupling matrix element between these two states 
must be an odd function of the in-plane momentum k. From these symmetry considerations, we 
deduce the general form of the effective Hamiltonian for the El and the HI states, expressed 
in the basis of \E1, mj = 1/2), \Hl,mj = 3/2) and \El,mj = -1/2), \Hl,mj = -3/2): 

H eff (k x , ky) = H ^_ k) ) , H = e(k) + d t {k)°i (2) 

where ai are the Pauli matrices. The form of H*(—k) in the lower block is determined from time 
reversal symmetry and H*(—k) is unitarily equivalent to H*{k) for this system(see Supporting 
Online Material). If inversion symmetry and axial symmetry around the growth axis are not 
broken then the inter-block matrix elements vanish, as presented. 

We see that, to the lowest order in k, the Hamiltonian matrix decomposes into 2x2 blocks. 
From the symmetry arguments given above, we deduce that d 3 (k) is an even function of k, 
while di(k) and d 2 (k) are odd functions of k. Therefore, we can generally expand them in the 
following form: 

d\ + id 2 = A(k x + iky) = Ak + 
d 3 = M- B(k 2 x + k 2 y ),e k = C- D{kl + fcj). (3) 

The Hamiltonian in the 2x2 subspace therefore takes the form of the 2 + 1 dimensional 
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Dirac Hamiltonian, plus an e(k) term which drops out in the quantum Hall response. The most 
important quantity is the mass, or gap parameter M, which is the energy difference between 
the El and HI levels at the T point. The overall constant C sets the zero of energy to be 
the top of the valence band of bulk HgTe. In a quantum well geometry, the band inversion in 
HgTe necessarily leads to a level crossing at some critical thickness d c of the HgTe layer. For 
thickness d < d c , i.e. for a thin HgTe layer, the quantum well is in the "normal" regime, where 
the CdTe is predominant and hence the band energies at the T point satisfy E(Tq) > E(Tg). 
For d > d c the HgTe layer is thick and the well is in the inverted regime where HgTe dominates 
and E(T 6 ) < E(T 8 ). As we vary the thickness of the well, the El and HI bands must therefore 
cross at some d c , and the gap parameter M changes sign between the two sides of the transition 
(Fig. 2). Detailed calculations show that, close to transition point, the El and HI band, both 
doubly degenerate in their spin quantum number, are far away in energy from any other bands 
(|2i]) . hence making an effective Hamiltonian description possible. In fact, the form of the 
effective Dirac Hamiltonian and the sign change of M at d = d c for the HgTe/CdTe quantum 
wells deduced above from general arguments is already completely sufficient to conclude the 
existence of the QSH state in this system. For the sake of completeness, we also provide the 
microscopic derivation directly from the Kane model using realistic material parameters (see 
SOM). 

Fig. 2 shows the energies of both the El and HI bands at hi = as a function of quan- 
tum well thickness d obtained from our analytical solutions. At d — d c ~ 64A these bands 
cross. Our analytic results are in excellent qualitative and quantitative agreement with previ- 
ous numerical calculations for the band structure of Hgi_-j.Cd-j.Te/HgTe/Hg! ^Cd^Te quantum 
wells (12711201) . We also observe that for quantum wells of thickness 40A < d < 70A, close to 
d c , the El± and Hl± bands are separated from all other bands by more than 30 meV (f2~71). 

Let us now define an ordered set of four 6-component basis vectors ipx 4 = (\E1, + > 
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, \H1, + >, \E1, — >, \H1, — >) and obtain the Hamiltonian at non-zero in-plane momentum 
in perturbation theory. We can write the effective 4x4 Hamiltonian for the El±, Hl± bands 
as: 

H?f f (k x , k y ) = [°° dz< vl)j\H{k x , ky, -id z )\iPi > . (4) 

J — oo 

The form of the effective Hamiltonian is severely constrained by symmetry with respect to z. 
Each band has a definite z symmetry or antisymmetry and vanishing matrix elements between 
them can be easily identified. For example, H^-J = 77g-^+ J^oo ^z(Tq, +~(z)\Tg, — \{z)) 
vanishes because |r 6 , +\){z) is even in z whereas |r 8 , —\){z) is odd. The procedure yields 
exactly the form of the effective Hamiltonian © as we anticipated from the general symmetry 
arguments, with the coupling functions taking exactly the form of ©. The dispersion relations 
(see SOM) have been checked to be in agreement with prior numerical results (|2ill20[) . We note 
that for k E [0, O.OlA -1 ] the dispersion relation is dominated by the Dirac linear terms. The 
numerical values for the coefficients depend on the thickness, and for values at d = 58A and 
d = 70A see SOM. 

Having presented the realistic k ■ p calculation starting from the microscopic 6-band Kane 
model, we now introduce a simplified tight binding model for the El and the HI states based 
on their symmetry properties. We consider a square lattice with four states per unit cell. The 
El states are described by the s-orbital states ^1,3 = \s,a = ±1/2), and the HI states are 
described by the spin-orbit coupled p-orbital states ^2,4 = i^bz ± Wy, a = ±1/2)- Here 
a denotes the electron spin. Nearest neighbor coupling between these states gives the tight- 
binding Hamiltonian of the form of Eq. [2} with the matrix elements given by 

d\ + id2 = A(sm(k x ) + i sin(fcj,)) 
d 3 = -25(2 - |§ - cos(k x ) - cos(ky)) 

e k = C -2D{2-cos{k x )-cos{k y )). (5) 
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The tight-binding lattice model simply reduces to the continuum model Eq. [2] when expanded 
around the T point. The tight-binding calculation serves dual purposes. For readers uninitiated 
in the Kane model and k ■ p theory, this gives a simple and intuitive derivation of our effective 
Hamiltonian that captures all the essential symmetries and topology. On the other hand, it also 
introduces a short-distance cut-off so that the topological quantities can be well-defined. 

Within each 2x2 sub-block, the Hamiltonian is of the general form studied in Ref. ©, in 
the context of the quantum anomalous Hall effect, where the Hall conductance is given by: 

a xy = -g-^ J J dk x dk y d ■ <9 x d x <9 y d (6) 

in units of e 2 /h, where d denotes the unit di(k) vector introduced in the Hamiltonian Eq. [2l 
When integrated over the full Brillouin Zone, a xy is quantized to take integer values which 
measures the Skyrmion number, or the number of times the unit d winds around the unit sphere 
over the Brillouin Zone torus. The topological structure can be best visualized by plotting d as 
a function of k. In a Skyrmion with a unit of topological charge, the d vector points to the north 
(or the south) pole at the origin, to the south (or the north) pole at the zone boundary, and winds 
around the equatorial plane in the middle region. 

Substituting the continuum expression for the di(k) vector as given in Eq. [3l and cutting 
off the integral at some finite point in momentum space, one obtains a xy = \sign(M), which 
is a well-known result in field theory (1221 . In the continuum model, the d vector takes the 
configuration of a meron, or half of a Skyrmion, where it points to the north (or the south) pole 
at the origin, and winds around the equator at the boundary. As the meron is half of a Skyrmion, 
the integral Eq. [6] gives ±|. The meron configuration of the di(k) is depicted in Fig. 2. In a 
non-interacting system, half-integral Hall conductance is not possible, which means that other 
points from the Brillouin Zone must either cancel or add to this contribution so that the total Hall 
conductance becomes an integer. The fermion-doubled partner (|23)) of our low-energy fermion 
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near the T-point lies in the higher energy spectrum of the lattice and contributes to the total a xy . 
Therefore, our effective Hamiltonian near the Y point can not give a precise determination of the 
Hall conductance for the whole system. However, as one changes the quantum well thickness 
d across d c , M changes sign, and the gap closes at the T point leading to a vanishing d^{k = 0) 
vector at the transition point d = d c . The sign change of M leads to a well-defined change of 
the Hall conductance Aa xy = 1 across the transition. As the di(k) vector is regular at the other 
parts of the Brillouin Zone, they can not lead to any discontinuous changes across the transition 
point at d = d c . So far, we have only discussed one 2x2 block of the effective Hamiltonian H . 
General time reversal symmetry dictates that a xy (H) = —a xy (H*), therefore, the total charge 
Hall conductance vanishes, while the spin Hall conductance, given by the difference between 
the two blocks, is finite, and given by AcjW = 2. From the general relationship between the 
quantized Hall conductance and the number of edge states (|2?|) . we conclude that the two sides 
of the phase transition at d = d c must differ in the number of pairs of helical edge states by 
one, thus concluding our proof that one side of the transition must be Z 2 odd, and topologically 
distinct from a fully gapped conventional insulator. 

It is desirable to establish which side of the transition is indeed topologically non-trivial. For 
this purpose, we return to the tight-binding model Eq. [51 The Hall conductance of this model 
has been calculated (24) in the context of the quantum anomalous Hall effect, and previously 
in the context of lattice fermion simulation (|25l) . Besides the Y point, which becomes gapless 
at M/2B = 0, there are three other high symmetry points in the Brillouin Zone. The (0, 7r) 
and the (%, 0) points become gapless at M/2B = 2, while the (ir, it) point becomes gapless at 
M/2B = 4. Therefore, at M/2B = 0, there is only one gapless Dirac point per 2x2 block. 
This behavior is qualitatively different from the Haldane model of graphene (1261) . which has 
two gapless Dirac points in the Brillouin Zone. For M/2B < 0, a xy = 0, while a xy = 1 for 
< M/2B < 2. As this condition is satisfied in the inverted gap regime where M/2B = 
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2.02 x 1CT 4 at 70A (see SOM), and not in the normal regime (where M/2B < 0), we believe 
that the inverted case is the topologically non-trivial regime supporting a QSH state. 

We now discuss the experimental detection of the QSH state. A series of purely electrical 
measurements can be used to detect the basic signature of the QSH state. By sweeping the 
gate voltage, one can measure the two terminal conductance Glr from the p-doped to bulk- 
insulating to n-doped regime(Fig. 3). In the bulk insulating regime, Glr should vanish at low 
temperatures for a normal insulator at d < d c , while Glr should approach a value close to 
2e 2 /h for d > d c . Strikingly, in a six terminal measurement, the QSH state would exhibit 
vanishing electric voltage drop between the terminals and fi 2 and between /x 3 and /i 4 , in 
the zero temperature limit and in the presence of a finite electric current between the L and R 
terminals. In other words, longitudinal resistance should vanish in the zero temperature limit 
with a power law dependence, over distances larger than the mean free path. Because of the 
absence of backscattering, and before spontaneous breaking of time reversal sets in, the helical 
edge currents flow without dissipation, and the voltage drop occurs only at the drain side of 
the contact (|77|) . The vanishing of the longitudinal resistance is one of the most remarkable 
manifestations of the QSH state. Finally, a spin filtered measurement can be used to determine 
the spin-Hall conductance ajgj . Numerical calculations (UH) show that it should take a value 
close tofjg = 2^. 

Constant experimental progress on HgTe over the past two decades makes the experimental 
realization of our proposal possible. The mobility of the HgTe/CdTe quantum wells has reached 
|i~6x 10 5 cm 2 /(Vs) (|27|) . Experiments have already confirmed the different characters of 
the upper band below (El) and above (HI) the critical thickness d c (12QII2ffl) . The experimental 
results are in excellent agreement with band- structure calculations based on the k -p theory. Our 
proposed two terminal and six terminal electric measurements can be carried out on existing 
samples without radical modification, with samples of d < d c ~ 64A and d > d c ~ 64A 



9 



yielding contrasting results. Following this detailed proposal, we believe that the experimental 
detection of the QSH state in HgTe/CdTe quantum wells is possible. 
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1 Figures 
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Figure 1 : (A) Bulk energy bands of HgTe and CdTe near the T point. (B) The CdTe/HgTe/CdTe 
quantum well in the normal regime El > HI with d < d c and in the inverted regime HI > El 
with d > d c . In this, and all subsequent figures r 8 /Hl (T 6 /E1) symmetry is correlated with the 
color red (blue). 
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Figure 2: (A) Energy (eV) of El (blue) and HI (red) bands at k\\ = vs. quantum-well 
thickness d (A). (B) Energy dispersion relations E(k x ,k y ) of the El, HI subbands at 40A, 
63. 5 A and 70 A from left to right. Colored shading indicates the symmetry type of band at that 
A; -point. Places where the cones are more red (blue) indicates that the dominant states are HI 
(El) states at that point. Purple shading is a region where the states are more evenly mixed. For 
40A the lower (upper) band is dominantly H1(E1). At 63. 5A the bands are evenly mixed near 
the band crossing and retain their d < d c behavior moving further out in A;-space. At d — 70A 
the regions near hi = have flipped their character but eventually revert back to the d < d c 
further out in A;-space. Only this dispersion shows the meron structure (red and blue in the same 
band). (C) Schematic meron configurations representing the di(k) vector near the T point. The 
shading of the merons has the same meaning a^the dispersion relations above. The change in 
meron number across the transition is exactly equal to 1, leading to a quantum jump of the spin 
Hall conductance Aa^j = 2e 2 /h. We measure all Hall conductances in electrical units. All of 
these plots are for Hgo.32Cdo.6sTe/HgTe quantum wells. 
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Figure 3 : (A) Experimental setup on a six terminal Hall bar showing pairs of edge states with 
spin up (down) states green (purple). (B)A two-terminal measurement on a Hall bar would 
give Glr close to 2e 2 /h contact conductance on the QSH side of the transition and zero on 
the insulating side. In a six-terminal measurement, the longitudinal voltage drops fi 2 — A*i and 
A*4 — A*3 vanish on the QSH side with a power law as the zero temperature limit is approached. 
The spin-Hall conductance a^j has a plateau with the value close to 2^-. 
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2 Supporting Online Material 

We show that our effective Hamiltonian can be derived perturbatively with k ■ P theory and is a 
quantitatively accurate description of the band structure for the El and the HI subband states. 
We start from the 6-band bulk Kane model which incorporates the T 6 and T 8 bands but neglects 
the split-off T 7 band( the contribution of the split-off band to the El, HI energies is less than 
5%(52)): 

H(k)=( EJ ^ + H < F/ T2X * ) (V) 

where E c is the conduction band offset energy, E v is the valence band offset energy, H C1 H v are 
the conduction and valence (Luttinger) band Hamiltonians, while T{k) is the interaction matrix 
between the conduction and valence bands: 







H r = ( ^ , 2fc2 ; Tt 







2m* 



! pi. 

— Ph .f^Ph 

o ±Pk„ ) 



^ = -£^ + H)k 2 + &M*-s) 2 ( 8 ) 



where m* is the effective electron mass in the conduction band, k± = k x ±ik y ,P = —^(s\p x \X) 
is the Kane matrix element between the s and p bands with m the bare electron mass, S is the 
spin-3/2 operator whose representation are the 4x4 spin matrices, and 71,72 are the effec- 
tive Luttinger parameters in the valence band. The ordered basis for this form of the Kane 
model is (|r 6 , +1/2 >, |T 6 , —1/2 >, |T 8 , +3/2 >,|r 8 ,+l/2 >, |T 8 , —1/2 >, |T 8 , —3/2 >). 
Although due to space constraints the above Hamiltonian is written in the spherical approxi- 
mation, our calculations include the anisotropy effects generated by a third Luttinger parameter 
73 / 72- The quantum well growth direction is along z with Hgx^Cd^Te for z < —d/2, HgTe 
for —d/2 < z < d/2 and Hg^^Cd^Te for z > d/2. Our problem reduces to solving, in the 
presence of continuous boundary conditions, the Hamiltonian Eq. [7] in each of the 3 regions 
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of the quantum well. The material parameters E c , E v , 71,2,3, and m*, are discontinuous at the 
boundaries of the 3 regions, taking the Hgi^Cd^Te/HgTe/Hgi ^Cd^Te values (see Ref. (SI) 
for numerical values), so we have a set of 3 multi-component eigenvalue equations coupled by 
the boundary conditions. The in-plane momentum is a good quantum number and the solution 
in each region takes the general form: 



7~t(k)ip(k x , ky, z) — l~C(k x ,ky, id z )ip{k X) ky ) z)\ 
ip(k x ,k y ,z) = e^+^^iz) 



(9) 



where ^(z) is the envelope function spinor in the six component basis introduced earlier. 

We solve the Hamiltonian analytically by first solving for the eigenstates at zero in-plane 
momentum, and then perturbatively finding the form of the Hamiltonian for finite in-plane k: 

H(k X) k y , —id z ) = H(0, 0, — id z ) + 5H(k x , k y , —id z ). At k x = k y = we have the following 
Hamiltonian: 



H(0,0,-id z ) 
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U + V J 
(10) 



where T = E c {z) + {-d z A{z)d z ), U = E v {z) - {-d zll {z)d z ), V = 2{-d zl2 {z)d z ). These 
parameters are treated as step functions in the ^-direction with an abrupt change from the barrier 
region to the well region. 

A general state in the envelope function approximation can be written in the following form: 



ty(k x ,k y ,z) = e 



i(k x x+k y y) 



fs(z) 

\Mz)j 



(11) 
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At k x = k y = the / 3 and / 6 components decouple and form the spin up and down (±) states 
of the HI subband. The fi, fn ft, h components combine together to form the spin up and 
down (±) states of the El and LI subbands. The linear-in-/c z operator ^Pk z , in H(0, 0, —id z ) 
forces the |T 6 , ±§)(z) and |r 8 , ±.\)(z) components of the El band to have different reflection 
symmetry under z <-> —z. The |r 6 ) band is symmetric in z (exponentially decaying in CdTe 
and a cosh(z) dependence in HgTe) while the |T 8 , mj = ±1/2) band is antisymmetric in z 
(exponentially decaying in CdTe and a sinh(^) dependence in HgTe). The opposite choice of 
symmetry under z — > —z reflection leads to the LI band. However, this band is far away in 
energy from both the El and the HI bands, does not cross either of them in the region of 
interest(52), and we hence discard it. 

For the El band we take the ansatz, already knowing it must be an interface state(52), to be: 
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(12) 



If we act on this ansatz with the Hamiltonian we decouple the 6 x 6 matrix into two, coupled, 
one-dimensional Schrodinger equations: 



Th{z) + ^-P( z )(-id z )U{z) = Eh{z) 
] -P{z){-id z )h{z) + (U - V)U{z) = Ef 4 (z) 



(13) 
(14) 



where P, T, U, V are given above. Using the restrictions from the Hamiltonian, the continuity 
of each wavef unction component at the boundaries, and the continuity of the probability cur- 
rent across the boundary, we derive the following set of equations that determine a and 5 as a 

function of E: 

E (Cd) _ A {Cd) a 2 {E) _ E _ \[\*a{E) 



2P 

3 i 



a{E) 



E (c d) + (7 po + 2l ic d))a2{E) _ E 
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E^-A^8 2 (E)-E = Vlfgg 

|f 5(E) ~ ' ~ + + 2 7 f 9) )5 2 (^) - E 

where a parameter X( Cd ) means the value of that parameter in the CdTe barrier material and 
X( H s) i s the value in the HgTe well material. Once we have a(E) and 8(E) we can use them to 
determine E through the following equation derived from the boundary conditions: 

E (cd) _ A {cd) a 2^ -E _ _ ^ 8(E)d ^ ^ E^) _ A iH ^8 2 (E) - (n) 



a(E) V 2 / V 6 ( E ) 

These rational transcendental equations are solved numerically to obtain the energy of the El 
subband at k x = k y = 0. 

We can follow a similar procedure to derive the energy of the HI subband. The heavy 
hole subband (at k x = k y = 0) completely decouples from the other bands and we have the 
one-dimensional, one-component Hamiltonian: 

Hf 3 (z) = E v (z) - ( 7l (z) - 2 l2 (z))(-d 2 z )f 3 (z) = Ef 3 (z). (18) 

We have the wavefunction in three regions 

\ ( C 3 e? z \ 

(19) 





( C*eP> \ 




V 3 cos(kz) 


\*in(z)J 


V C 3 e-P* J 



where f3 2 (E) = and k 2 (E) = { ^ g) ^ Hg) . We can pick C 3 = 1 which gives us the 

relation 

Vs = . (20) 
cos(K(E)a/2) 

from the boundary condition at z = —d/2. Finally, we need to normalize the wavefunction to 
get the coefficients. The energy of this state is determined by considering the conservation of 
probability current across the boundary. The following equation is solved for the energy: 

' cot(«(f?)d/2). (21) 



(^-24 Cd) )P(E) (^-2^)n(E) 
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We repeat both of these processes on a state with only f<z(z) and f 5 (z) non-zero and on a 
state with only fe(z) non-zero and to get the E\~ and Hl~ bands respectively. We have the 
forms of these states at k x = k y = and can use k ■ P perturbation theory to derive a two 
dimensional Hamiltonian near the T point in /c-space. 

3 Perturbation Theory and Effective Hamiltonian 

Define an ordered set of basis vectors (\E1, + >, \H\, + >, \E1, — >, \H\, — >). We can 
write the effective Hamiltonian as: 

/oo 
dz < ^j\HQ X Q{k x , k y , -id z )\ipi > (22) 
-oo 

where ipi is the i-th element of the basis set given above which will give a 4 x 4 effective 
Hamiltonian. The integrals must be split into the three regions defined above, and the parameters 
from each material must be accounted for in the Hamiltonian. This Hamiltonian depends on 
the quantum-well width d, and once d is specified we can numerically calculate the matrix- 
elements. It is important to note that f\(z), f2,(z), fs(z), fe(z) are symmetric with respect to z 
and f 4 (z), f*,(z) are antisymmetric in z(S2) which is a useful simplification in performing the 
integrals. An example of one integral is of the form: 

Jdzf*(z) { l3 (z),-td z }f A (z) = 
J dz 1 - (2 l3 (z)f* 3 (z)d z f 4 (z) + K(z)f A (z)dM*)) ■ (23) 
The functional form of 73(2) is 

7s (z) = 7 f d) (6(-d/2 -z) + 6(z - d/2)) + 7 f fl) (0(* + d/2) - 6{z - d/2)). (24) 

The ^-derivative acting on this function produces 5-function terms that contribute a term pro- 
portional to (73^ — li Hg ^) to the integral which vanish when these material parameters are 
equal. The integral then has to be evaluated numerically. 
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d(A) 


A(eV) 


B(eV) 


C(eV) 


D(eV) 


M(eV) 


58 


-3.62 


-18.0 


-0.0180 


-0.594 


0.00922 


70 


-3.42 


-16.9 


-0.0263 


0.514 


-0.00686 



Table 1: Parameters for Hgo.32Cdo.6sTe/HgTe quantum wells. 



After calculating the matrix-elements we are left with an effective Hamiltonian parameter- 
ized in the following way: 

\ 



7~t{k x , ky) 



fe k + M(k) Ah. 

Ak + e k -M{k) 

e k + M{k) -Ak 

V -Ak. e k -M(k)J 



(25) 



where e k = C - D(k 2 x + fcj), M(k) = M - B(k 2 x + k 2 ), k± = k x ± iky, and A, B, C, D, M 
depend on the specified quantum-well width. For values of these parameters at d = AO A and 
d = 70 A see Table 1. This Hamiltonian is block diagonal and can be written in the form 

H(k) 



7i,[k x , ky) — 



(26) 



H*(-k) 

where H(k) = e k I 2X 2+d a (k)a a ,withd 1 = Ak x ,d 2 = Ak y ,andd 3 = M(k) = M-B(k 2 x +k 2 y ). 
Finally, we define a unitary transformation: 



U 



hx2 
-o 1 



(27) 



and take WH 4x4 U which reverses the sign of the linear k terms in the lower block and puts it 
H into the form 



7~t(k x , ky) 



H(k) 



(28) 



H*(k) 

U does not affect the ^-direction of spin and simply rotates the x and y axes in the lower block 
of the Hamiltonian by ir. The energy dispersions for these bands are given in Fig. 1 of the 
supporting online material for several values of d. 
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Figure 4: Dispersion relations for the El and HI subbands for(A) d = 40A (B) d = 63. 5A (C) 
d = 70A. 
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